Multichannel Generalization of Kitaev's Majorana End States and a Practical Route 

to Realize Them in Thin Films 
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The ends of one-dimensional p+ip superconductors have long been predicted to possess localized 
Majorana fermion modes[l|]. We show that Majorana end states are robust beyond the strict ID 
single-channel limit so long as the sample width is not much larger than the superconducting co- 
herence length, and exist when an odd number of transverse quantization channels are occupied. 
Consequently the system undergoes a sequence of topological phase transitions driven by changing 
the chemical potential. These observations make it feasible to implement quasi- ID p+ip super- 
conductors in metallic thin-film microstructures, which offer 3-4 orders of magnitude larger energy 
scales than semiconductor-based schemes. Some promising candidate materials are described. 

PACS numbers: 71.10.Pm, 74.20.Rp, 74.78.-w, 03.67.Lx 



Since the pioneering work by Kitaev[l], there has been 
growing interest in the condensed matter community to 
find a way to realize Majorana fermions in the laboratory. 
Often described as "half a fermion" , Majorana fermions 
are expected to exhibit non-Abelian braiding statistics 
[2|, |3[ and have been proposed as a basis for topologi- 
cal quantum computers which would be protected from 
decoherencejjj, |^. While numerous schemes for realiz- 
ing Majorana fermions have been proposed |5l4l5|, up 
to now, all face significant technical challenges to imple- 
ment. The early proposal to use the v = 5/2 fractional 
quantum Hall state |5| has not yet received definite con- 
firmation despite significant experimental efforts. The 
proposal by Fu and Kane J6| to use proximity induced 
superconductivity in the surface state of topological in- 
sulators opened up a promising new avenue but awaits 
further development in topological insulator materials. 

Recently, a number of proposals based on more tra- 
ditional materials have appeared. The basic idea is to 
use the Rashba-type spin-orbit coupling in combination 
with conventional 5-wave superconductivity to produce 
an effective Px ± ipy 2D superconductor. Magnetiza- 
tion is then introduced to remove one of the two com- 
ponents, leaving an effective px + ipy superconductor J7l- 
Il0| which is known to carry a single Majorana fermion 
in its vortex cores [2|, |3|. The Rashba coupling may be 
realized in semiconducting quantum wells [sl, l9[, non- 
centrosymmetric bulk superconductors, or on the sur- 
faces and interfaces of strong spin-orbit coupled con- 
ventional superconductors [71 |lO|. In one scheme a sin- 
gle component Px + ipy superconductor would be real- 
ized by using the spin splitting of a ferromagnetic metal 
thin filmjlO|. Unfortunately, this requires control of the 
film thickness and interfaces to atomic precision in order 
to achieve an odd number of modes. A more sophisti- 
cated scheme uses the fact that out-of-plane magnetiza- 
tion opens up a gap A^ at the Dirac crossing of a Rashba 
band. As long as Am is greater than the superconducting 
gap A, an effective px + ipy superconductor results [7H9|. 



However, the electron density in such structures is ex- 
tremely low [9], and exquisite purity will be needed to 
avoid carrier localization. 

In addition to the above-mentioned challenges facing 
the construction of an effective px + ipy superconduc- 
tor, realizing Majorana fermions as bound states in vor- 
tex cores presents the additional difficulty of controllably 
creating and manipulating multiple superconducting vor- 
tices. This has led to renewed interest in the original 
Kitaev idea |12l - ll4| , where Majorana particles are shown 
to be localized at the ends of a one-dimensional Px + ipy 
superconducting wire. In an interesting recent paper, 
this setup has been examined in great detail, including 
explicit plans to manipulate, braid, and measure these 
Majorana end states [15|. 

All previous discussions of Majorana end states, 
with one exception [l7|, assume that an effective one- 
dimensional band can be created by occupying only the 
lowest transverse sub-band of a quantum wire. This 
again presents formidable technical challenges. In this 
paper we move beyond the strict one-dimensional limit, 
and explore further the notion of Majorana end states in 
quasi-one dimensional p^^ + zp^/ superconductors. We con- 
sider samples of finite width W, and allow the occupation 
of multiple transverse sub-bands. Our key finding is that 
the Majorana end states are surprisingly robust. They 
exist when an odd number of transverse channels are oc- 
cupied, and are protected from low lying excitations as 
long as the sample width W does not substantially exceed 
the superconducting coherence length ^o = ttvf/A. Thus 
a series of quantum phase transitions between topological 
and conventional states occurs as a function of chemical 
potential. Due to the condition W < ^o^ the Majorana 
end states are stable over a range in chemical potential 
which exceeds A. We also find that these zero modes 
survive disorder, as long as the mean free path £ exceeds 
W. 

These results suggest that it is feasible to create 
Majorana states in thin films. To this end it is necessary 



to create a 2D Px^iPy superconductor. In the second half 
of this paper we propose a new scheme which combines el- 
ements of earlier proposals [7|, |8|, llO| and uses the surface 
states of a strong spin-orbit coupled metal with super- 
conductivity and magnetization induced by proximity ef- 
fect. The advantage over semiconductor heterostructure 
schemes is that the energy scale of the Rashba splitting 
and the electron density are much higher. Some promis- 
ing candidate materials are described. 

Majorana end states in quasi- one- dimensional wires - 
We consider a rectangular px + ipy superconductor of 
length L ^ ^0, and width W . For sufficiently large W, 
the system is fully two-dimensional and the low lying 
{e < A) states are chiral edge modes with a gap set by the 
finite circumference of the sample. This gap is required 
by the fact that a single zero-energy Majorana state can- 
not exist in isolation, but can only occur in pairs. The 
question is: how are these chiral edge states connected 
with the zero energy Majorana end states as the width 
decreases? 

To probe this crossover from 2D to ID, we consider a 
square-lattice tight-binding model, H = Ht -\- i^p-BCS, of 
a single species of electrons with px + ipy BCS pairing: 



Ht = '^-t (clcj + h.c. j - Y^ ^c]cj 

ifp.Bcs = E ^ {-^4+^4 + 4+ A) + h-c (1) 



where cj creates an electron on site j, t is the hopping 
amplitude, /i is the chemical potential, A is the p-wave 
pairing amplitude, and we work in units where the lattice 
spacing is unity. 

Fig. [U shows the lowest energy state wave-function, 
obtained by numerically diagonalizing Eq. ([!]), for a se- 
quence of samples with increasing widths. For W < ^o, 
there are two well isolated Majorana end states local- 
ized at opposite ends. As W is increased to ^ ^o, these 
Majorana end states begin to spread along the edges of 
the sample with extent ~ e^^^^ . For fixed L, these wave- 
functions eventually circumnavigate the sample to match 
up with the 2D edge mode picture for W ^ ^q. The re- 
sults of these simulations can be understood as follows: 
in the 2D limit, the edge modes are localized near the 
sample boundary and decay into the bulk with charac- 
teristic length ^0- As VK is decreased and approaches ^o, 
the tails of the edge mode wave-functions on opposite 
edges overlap and mix. The mixing of these two counter- 
propagating edge modes forms a gap of order e~^^^^ 
along the mid-section of the sample, and squeezes the 
lowest energy state towards opposite ends of the sample. 
For sufficiently long samples {L ^ e^^^°) one finds spa- 
tially isolated zero-energy Majorana end states residing 
in opposite ends. 

As in shown Fig. [2li, the Majorana end modes are pro- 
tected from low-energy excitations by a spectral gap that 
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FIG. 1. (Color online) |^|^ for zero-mode wave- functions at 
fixed coherence length: ^o ^ t/A = 10, chemical potential 
/i = — 2t, fixed sample length: L = 200, and various sample 
widths: W = 10, 30, 50 and 100 (from left to right). These 
widths span the range from W < <Jo, W ^ ^o, to W ^ ^o, 
and depict the crossover from the ID regime with Majorana 
zero- modes localized at each end, to the 2D limit with chiral 
edge mode. 



scales as ~ Ae~^^^^. Therefore, while a sample of any 
W can in principle support Majorana end-states for suf- 
ficiently large L ^ e^Z-^o ^ [j^ practice we have the restric- 
tion W ^ ^0 to avoid an exponentially small excitation 
gap. For W < ^0 and L :^ (_o^ this gap is a sizeable frac- 
tion of the bulk superconducting gap. A, and is largely 
insensitive to both L and W. 

Topological phase transitions - In the quasi- ID regime, 
where W ^ ^o, the Majorana end states discussed above 
occur only when an odd number of transverse sub-bands 
is occupied. We next examine the effects of changing 
chemical potential /i. The top panel of Fig|2] shows the 
low lying excitation energies as a function of fi for an 
L = 100, W = 10 lattice with t/A = 10-^0- The 
higher energy states shown in blue reside in the bulk 
and form two branches: excitations across the pairing 
gap and excitations to the next higher sub-band which 
respectively increase/decrease in energy with /i. Due to 
the p-wave pairing symmetry, the pairing gap vanishes 
when fi coincides with a sub-band bottom (dashed hues). 

We now focus on low energy states lying within the 
bulk gap. Starting from the lowest transverse sub-band 
and increasing /i, one finds an alternating sequence of 
transitions between states with a non-degenerate zero- 
mode (shown in red) and doubly degenerate gapped 
states as described above. The transition points coin- 
cide with /i passing through the bottom of a transverse 
sub-band, and are accompanied by a closing of the bulk 
gap (as is required for any topological phase transition). 

Intuitively one can imagine that each sub-band con- 
tributes a Majorana mode at each end of the sample. An 
even number of Majorana fermions localized at a given 
end will mix and pair into full electron states at non-zero 
energy. For an odd number, however, there is always one 
remaining unpaired Majorana mode per end. Thus, as a 
function of chemical potential /i the system undergoes a 
series of topological phase transitions between topolog- 
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FIG. 2. (Color online) a) Excitation spectrum as a function of 
chemical potential, /i, for 100 x 10 site lattice strip of p^ -\-ipy 
superconductor with t/A = 10 ^ .Jo- Blue lines are bulk 
excitations. Dashed lines denote the transverse confinement 
band-bottoms, which coincide with topological phase transi- 
tions between phases with Majorana zero-modes (shown in 
red), to topologically trivial gapped phases, b) and c) the 
effect of random on-site bulk disorder, normally distributed 
with width Wn/t — 50%, and 75% (left to right respectively). 
The Majorana zero-mode survives moderate disorder, but is 
destroyed for sufficiently strong disorder. Insert d) shows the 
exponential sensitivity of the excitation gap Aex to W for 
clean samples with L — 1000, t/A = 5 '^ ^o and fi near — 2t 
(fine tuned to support a zero- mode). 



ical states with spatially isolated Majorana zero-modes 
to topologically trivial gapped states with no Majorana 
modes. We have verified this picture by checking that the 
non-degenerate (red) states have energy that is exponen- 
tially small in L, and that their wavefunctions can be 
written as ?/^ = a + z6, where a and h are real fermions lo- 
calized in opposite ends. On the other hand, the gapped 
states are doubly degenerate due to spatial parity symme- 
try {x^y) -^ —{x^y) and are complex fermions localized 
in the sample ends. 

Robustness against disorder - In realistic experimen- 
tal devices, some degree of random disorder is unavoid- 
able. To explore the effect of bulk disorder, we add a 
random onsite potential i^Disorder = ^j ^j^j^j to the 
tight-binding model ([1]), where {Vj} are normally dis- 
tributed with variance Wd- Fig. [2]3,c show a segment 
of the excitation spectrum as a function of /i for typi- 
cal disorder realizations at two disorder strengths. For 
weak to moderate disorder strength, the Majorana zero- 
mode survives, and the main effects of disorder are to 
narrow the spectral gap to the lowest excited state and to 
break the parity degeneracy of the excited states. How- 
ever, for sufficiently strong disorder, the excitation gap 
closes and the Majorana modes no longer exist for any 
ji. We have checked that the critical disorder strength 
depends on the relative size of the sample width and 
the mean- free-path £ = vp [7rW^N{/jj)) (where N{e) 
is the density of states, and vp is the Fermi velocity) 




FIG. 3. (Color online) Proposed metallic thin film based struc- 
ture for creating apx-\-ipy superconductor (a). A surface state 
with large Rashba splitting has been observed in the Ag(lll) 
surface (gray circles) covered by 1/3 monolayer of Bi (purple 
circles) [l^. Adjacent layers of superconducting Nb alloy and 
ferromagnetic insulating EuO induce single-species px + ipy 
pairing in the surface state. Rashba and ferromagnetic split 
band structure (b) and corresponding density of states, N(€) 
(c). N{e) has a Van-Hove singularity at the Rashba band- 
bottom, and exhibits a dip due to the ferromagnetic magne- 
tization gap Am = 2Vz. 



with the Majorana mode surviving for i ^ W. It is im- 
portant to note that, for moderate disorder strength, the 
existence of Majorana-modes and location of topological 
phase transitions are no longer precisely tied to the trans- 
verse quantization band-bottoms. This suggests that the 
existence of Majorana zero- modes and alternating even- 
odd phase transition structure do not rely crucially on 
sharp transverse quantization. 

Another potentially important source of disorder is 
spatially varying edge orientation. To assess its impact, 
we have also performed simulations for strips with me- 
andering piece- wise-linear edges. We find similar results 
to those of bulk-disorder, with Majorana end-states sur- 
viving moderate spatial variations in iy[16]. 

Thin-film realization - In order to realize the Ma- 
jorana end-states discussed above, one must first con- 
struct di Px -\- ipy superconductor. As proposed in [7|- 
9], an effective Px + ipy superconductor can be engi- 
neered by combining three ingredients: 1) a 2D sur- 
face band with strong Rashba spin-orbit coupling term 
^R = o^rz • (<j X k), 2) s-wave spin-singlet pairing, and 
3) a ferromagnetic splitting term, H^u = VzCFz^ to effec- 
tively remove one 'spin-species'. We have shown above 
that Majorana end states in quasi- ID Px + ipy super- 
conductors survive beyond the strict ID limit and are 
robust against disorder. These features make it feasi- 
ble to realize strong Rashba coupling in the above recipe 
using metallic thin films rather than semiconductor het- 
erostructures. 

A particularly promising example is the surface state 
which forms when 1/3 monolayer of high Z elements such 
as Bi is grown on Ag(lll) surface. This surface state is 
well isolated from bulk states and exhibits large Rashba 
coupling on the order of A^^ = 2aRkF — ^vffo?^ = 



1.85eV[l8|. This value exceeds typical numbers achiev- 
able in semiconductor structures by three to four orders 
of magnitude^. As discussed in [9|, large A^^ not only 
ensures that the induced s-wave pairing is effectively con- 
verted into px + Wy pairing, but also determines the 
maximum achievable carrier density in the topological 
superconducting phase. These considerations lead us to 
propose the structure shown in Fig. [3l 

In this scheme, electrons confined to the Bi/Ag surface 
form a purely 2D Rashba split surface band. Supercon- 
ductivity can be induced by proximity to a conventional 
superconductor such as Nb. Ferromagnetic splitting can 
be achieved by depositing a layer of ferromagnetic insu- 
lator, such as EuO, on top of the Bi/Ag surface layer. 
Experiments on ultra-thin films of EuO have demon- 
strated sizeable ferromagnetic exchange energies equiv- 
alent to Vz ^ 300meV[l9|. While the actual induced 
Zeeman splitting will depend sensitively on the trans- 
parency of the EuO/Bi-Ag interface, even a small frac- 
tion of this previously demonstrated Vz would be suffi- 
cient for our purposes. Planar tunneling could be used 
to characterize the surface Rashba and induced Zeeman 
energies. Fig. [3fc shows the density of states N{e) for a 
surface band with Rashba and Zeeman terms. N{e) has 
a square-root van Hove singularity at the Rashba band- 
bottom which has been observed directly in tunneling 
data[20|. N{e) also exhibits a dip due to the magnetiza- 
tion gap A^ = 2Vz. This dip is present only when the 
induced magnetization lies perpendicular to the surface, 
and could potentially be resolved by comparing tunnel- 
ing data for in-plane and perpendicular magnetization. 
Finally, in order to induce single-species Px + Wy pair- 
ing, ji must be tuned to lie within the Zeeman gap (see 
Fig. [3b). This tuning can be achieved, for example, by 
electrostatic gating. Local tunneling experiments can be 
used to unambiguously verify the presence of Major ana 
modes. These would show up as zero-energy mid-gap 
states, and at sufficiently low temperature would exhibit 
quantized conductance of G = 2e^ /h for resonant An- 
dreev reflection [21 1. 

We have also carried out simulations of finite strips 
of spin-full fermions in the presence of Rashba, Zeeman 
and s-wave pairing terms. The results agree with those 
shown above, and display the alternating appearance and 
disappearance of Majorana fermions with /i, as long as /i 
lies within the Zeeman gap. 

Conclusions and Discussion - We have shown that 
robust Majorana end states in ID px + ipy supercon- 
ductors survive beyond the strict ID limit, so long as 
the sample width W is not significantly larger than the 
superconducting coherence length ^o, opening the door 
for realizing Majorana fermions in metallic thin film 
microstructures with large Rashba spin-orbit coupling. 
Such materials offer orders of magnitude larger energy 
scales and carrier concentrations compared to analogous 
semiconductor-based proposals. The Majorana states 



can be controlled either by fine-tuning ji with electro- 
static gates, or by controlling domain walls in the Fer- 
romagnetic insulator. For example, an effective Px + ipy 
region can be created at a domain wall separating two 
regions of opposite in-plane magnetization. Our results 
raise the exciting possibility that Majorana fermions 
could be realistically implemented in microstructures de- 
fined by conventional lithography techniques, and could 
be efficiently manipulated and braided simply by electro- 
static gating or other means. 
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Appendix A: Simulations with a Meandering 
Edge: In addition to random bulk on-site disorder, one 
may also worry about the effects of spatially varying non- 
parallel edges. This type of edge disorder mixes trans- 
verse sub-bands by edge-scattering, and if sufficiently 
strong may close the bulk gap and destroy the Majo- 
rana end- states. However, we find that the Major ana 
end-states are robust against moderate spatial variations 
in edge width and orientation. 

To demonstrate the robustness of the multi-channel 
scheme to moderate amounts of edge- variance, we have 
conducted simulations of p -\- ip strips with meandering 
edge profiles. Random meandering edges were created by 
introducing a strong electrostatic confinement potential 
Konf(^,^), which is zero inside the strip and transitions 
to V^conf (^, y) = 4t outside the strip. The top and bot- 
tom edges were independently formed by dividing the 
length of the strip into uniformly random intervals of av- 
erage length Lmeander- The positiou of the top and bot- 
tom edges were chosen randomly at the endpoint of each 
of the aforementioned intervals with average ±W/2 and 
Gaussian spread aw • In between these end-points, the 
top and bottom edge profiles were linearly interpolated 
(in discrete steps). Finally, in order separate the effects 
of non-parallel edges from edge jagged- ness, Vconf{y) was 
smoothed over the length scale of one lattice spacing by 
applying a Gaussian filter along y. 

Representative simulations of samples with meander- 
ing edges are shown in Figures |4] and [H Fig. |4]a.-d. show 
colormaps of Konf for increasing edge variance: aw = 0, 
1, 1.5, and 2 lattice spacings respectively. These vari- 
ances represent fractional variances in the sample width 
of 0%, 6.7%, 10%, and 13.3% respectively. 

Fig. [5] shows the corresponding excitation spectra, 
as a function of chemical potential /i. In these figures, 
energies are measured with respect to typical straight- 
edge excitation gap scale of Aex^^ = O.IA shown in (a). 
Like bulk-disorder, the effect of edge disorder is to re- 
duce the bulk excitation gap, with larger variance in the 
edge geometry (either through increasing aw or decreas- 
ing i^meander) leading to Stronger reduction. 

Crucially, the Majorana zero-modes (highlighted in 
red) survive moderate amounts of edge- variance. Fur- 
thermore, the alternating sequence of topological phase 
transitions as a function of fi survives with meandering 
edges. However, as with bulk disorder, the locations of 
the phase transitions are no longer tied to the clean trans- 
verse sub-band bottoms. The amount of edge- variance 
that can be tolerated depends on the system geometry 
(W, L, and ^meander) and coherence length ^o- From sim- 
ulations of samples with various widths (but with fixed 
^o/W ~ 1), we find that typically, the Majorana end- 



states can survive up to ^ 15% variation in the edge 
width. 
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FIG. 4. (Color online) Meandering edge confinement poten- 
tial, Konf, spatial profiles for 300x30 strip with t = 30, A = 1. 
(a)-(d) show random samples with fixed meandering corre- 
lation length: Lmeander = 60 = 3VK, and increasing width 
variance: aw = 0, 1, 1.5, and 2 respectively. Corresponding 
excitation spectra shown in Fig. [5] below. 
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FIG. 5. (Color online) Corresponding excitation spectrum as 
a function of /i, for meandering edge samples shown in Fig. 
IHa)-(d) above. Excitation energies are measured with re- 
spect to the typical excitation energy Aex^^^ = O.IA for the 
straight-edge case (a). The effects of the smooth edge disor- 
dering are similar to those of bulk onsite disorder: increasing 
aw decreases the bulk excitation gap by mixing transverse 
sub-channels. However, for moderate aw the Majorana end- 
states survive and continute to appear and dissapear in an 
alternating fashion as /i is swept. 



